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ABSTRACT 

The leading long-distance quantum correction to the Newtonian potential for 
heavy spinless particles is computed in quantum gravity. The potential is obtained 
directly from the sum of all graviton exchange diagrams contributing to lowest non- 
trivial order to the scattering amplitude. The calculation correctly reproduces the 
leading classical relativistic post-Newtonian correction. The sign of the perturbative 
quantum correction would indicate that, in the absence of a cosmological constant, 
quantum effects lead to a slow increase of the gravitational coupling with distance. 
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1 Introduction 



It is generally assumed that a quantum theory of gravity cannot lead to testable 
predictions, due to a lack of perturbative renormalizability of the Einstein-Hilbert 
action [1, 2, 3, 4, 5]. Recently the interesting possibility has been raised [6] that 
low energy predictions of quantum gravity are not necessarily affected by the short 
distance details of an ultraviolet regulated theory of gravity [7], and can in fact be 
finite and calculable. As an application, the leading long distance quantum correc- 
tion to the static Newtonian potential was computed, resulting in a finite correction 
of order 0(Gh / c 3 r 3 ). When gravity is treated in this fashion as an effective low 
energy theory, the analogy with the use of effective field theories in treating the 
physics of soft pions comes to mind [8]. 

The existence of a universal long distance quantum correction to the Newtonian 
potential should be relevant for a wide class of gravity theories. It is known that 
the ultraviolet behavior of pure Einstein gravity can be improved by adding higher 
derivative contributions to the action [7]. In four dimensions these can be restricted 
to the form aR^R^ + /3i? 2 , where a and f3 are dimensionless coupling constants. 
The resulting classical and quantum corrections to gravity are expected to alter 
significantly the potential at short distances (comparable to the Planck length) but 
should not affect the behavior at large distances, which should largely be deter- 
mined by the structure of the Einstein-Hilbert action. Only the latter action will 
be therefore the subject of our present investigation. We should add that simpli- 
cial lattice regularizations of gravity also include in general higher derivative terms, 
and the same considerations should apply in this well, as long as the correct 

continuous gauge invariance properties of the continuum action are incorporated [9]. 

In the following we will compute the leading classical and quantum corrections to 
the static potential, by evaluating the complete set of diagrams which contribute to 
the scattering amplitude for heavy spinless particles in the low momentum transfer 
limit. From the resulting expression the effective static potential at large distances 
can then be read off easily, and will contain, as explained further below, both classical 
relativistic and quantum corrections. An important omission in our calculation 
will be the absence of a bare cosmological constant, which would complicate the 
perturbative treatment significantly due to the need to expand about a non-flat 
background. Our results for the static potential are such that they suggest a slow 
increase of gravitational interactions with distance due to the quantum correction. 
The answer will then be compared with two recent calculations of the same quantity 
[6, 10, 11], which include only a subset of the diagrams considered here. We will find 
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that our answer is qualitatively similar to the result of [6], but differs in sign from 
the result of [10], where a rather different method, based on world-line correlations, 
is used to estimate the potential. 

2 One Loop Amplitudes 

Before describing the calculation, it will be useful to first clarify our conventions 
and notation. We shall expand around the flat Minkowski space-time metric, with 
signature given by rj^ = diag(l, — 1, — 1, — 1). The Einstein-Hilbert action is then 
given by 

Se = + T^g I dx \/-9( x )R( x ), (2-i) 

with g(x) = det(g fll/ ) and R the scalar curvature. It is also assumed in the following 
that the bare cosmological constant is zero. The presence of a non- vanishing cosmo- 
logical constant introduces additional momentum independent vertices, and would 
make the perturbative calculation described below considerably more difficult. In 
particular the expansion around flat space is no longer justified in this case, and 
it has to be performed around a solution of Einstein's equations with a non-zero 
cosmological constant. 

The coupling of gravity to scalar particles of mass m is described by the action 

S m = l -j 'dx y/-g(x) [g^(x)d^(x)cU(x) ~ m 2 ^ 2 (x) ] (2.2) 

In the following we shall consider the interaction induced by graviton exchange 
between two heavy scalar particles of distinct mass m\ and m 2 . The effective inter- 
action in the static limit is then determined by evaluating the scattering amplitude 
between the two heavy particles, in the limit of small momentum transfer q 2 — > 0. 

Usually in perturbation theory the metric g llv {x) is expanded around the flat 
metric rj^ [4], by writing 

g^{x) = j}^ + k h^x) (2.3) 

with k, 2 = 32ttG. Here we shall instead follow the method of reference [12], and 
define the small fluctuation graviton field h fll/ (x) via 

g^{x)^-g{x) = V ^ + nh^(x) (2.4) 

One advantage of this expansion over the previous one is that it leads to considerably 
simpler Feynman rules, both for the graviton vertices and for the scalar-graviton 
vertices (as can be seen from the fact that precisely the above expression appears 
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in the kinetic term of the scalar field action). Once the action is expanded out in 
the graviton field h^ v {x), the space-time indices are then raised and lowered using 
the flat metric, and there is therefore no longer a need to distinguish between upper 
and lower indices. 

A gauge fixing term [13, 14] has to be introduced, and here it will be of the form 

^(d,^{x)g^)\ (2.5) 
as suggested in Ref. [12]. The bare graviton propagator is then given simply by 

D^ P a{p) = j— • 2.6 

2 p z + it 

For the present calculation one also needs expressions for the three-graviton and 
two ghost-graviton vertex. The relevant expressions are quite complicated and have 
already been given in Ref. [12], so they will not be reproduced here. We have 
performed a number of checks of the results of Ref. [12], some of which will be 
discussed below. Let us point out here that with the present definition for the 
gravitational field, there are no factors of l/(d — 2) for the graviton propagator in 
d dimensions; such factors appear instead in the expressions for the Feynman rules 
for the vertices. For the following calculations we shall also need the two scalar-one 
graviton vertex, which is given by 

IK i ^ 2 



2 yPinP2v + PivP2f> ~ m r)^j (2.7) 

where the P\ } p2 denote the four-momenta of the incoming and outgoing scalar field, 
respectively. In addition we need the two scalar-two graviton vertex, which is given 
by 

ik m 



2(d - 2) + ~ ~d~2 J ( 2-8 ) 

where one pair of indices (//, v) is associated with one graviton line, and the other 
pair (A, a) is associated with the other graviton line. These rules follow readily from 
the expansion of the gravitational action to order G 3 / 2 (k 3 ), and of the scalar field 
action to order G (k 2 ). 

To lowest order in G, the contribution to the potential can be computed from 
the single graviton exchange diagram. In momentum space the static contribution 
is given, as expected, by 

-Gm x m 2 % (2.9) 
q z 

where q is the momentum transfer (see also [15]). 
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Higher order corrections in G are computed by evaluating contributions to the 
interaction coming from the complete set of one-loop diagrams. One notices that 
the relevant length scale appearing with the Einstein-Hilbert action for pure gravity 



particle involves only the combination mc/h } the inverse Compton wavelength asso- 
ciated with the heavy sources. This is also clearly seen from the path integral phase 
contribution for a single particle, which is given by 



When one considers the lowest order contribution to the gravitational interaction 
due to single graviton exchange one obtains a contribution to the static gravita- 
tional potential proportional to {%/ c)(mc 2 /h) 2 (Gh/ c 3 ) = m 2 G. At order G 2 one 
finds contributions both of order {%/ c)(mc 2 /h) 2 (Gh/ c 3 ) 2 = m 2 fiG 2 / c 3 and of or- 
der (h/ c)(mc 2 /H) 3 (Gh/ c 3 ) 2 = m 3 G 2 /c. The first one represents a genuine quantum 
correction proportional to h, while in the second type of contribution the /j's have 
canceled, and the resulting correction represents a classical relativistic correction. 
The latter involves the Schwarzschild radius of the massive particle, 2Gm/c 2 . 

These considerations lead to the apparently paradoxical result that Feynman 
diagram perturbation theory is also expected to reproduce the classical relativistic 
corrections, which are independent of fi. Indeed it was shown by the authors of 
Ref. [16, 17] that classical relativistic corrections involve tree graphs connected to 
an arbitrarily high number of external classical sources. A detailed calculation of 
these classical relativistic corrections, using diagrammatic methods, was performed 
in Ref. [18]. There it was shown explicitly that the corrections of order G 2 correctly 
and completely reproduce the leading classical relativistic corrections appearing in 
the Einstein-Hoffmann- Infeld effective post-Newtonian Hamiltonian. 

Let us now return to the computation of the one loop amplitude. One needs 
to calculate all first-order corrections in G, which will include both the classical- 
relativistic 0(G 2 m 2 / c 2 ) and the quantum mechanical 0(hG 2 / c 3 ) corrections to the 
classical Newtonian potential energy discussed before. The relevant topologically 
distinct Feynman diagrams are shown in Figs. 1 and 2. 




On the other hand the action for the scalar 
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la) 



lb) 





lc) 



Id) 



Fig 1. Some one loop graviton exchange diagrams. 

The relevant amplitudes are computed in momentum function of the 

total momentum transfer squared q 2 . They are then evaluated using dimensional 
regularization in 4 — e dimensions, using the Feynman parametric representation for 
combining propagator denominators. The final answer then follows after performing 
the necessary momentum and parametric integrations. Due to the vast amount of 
algebraic manipulations involved in doing the index contractions, computer algebra 
was employed throughout the calculation in order to ensure the correctness of the 
results. For small q 2 the contributions arising from each diagram can then be sep- 
arated into two types of terms, one describing the classical relativistic correction 
proportional to and the other describing the leading quantum correction 

proportional to Xogq 2 . 

These in turn can then be expressed as corrections in coordinate space by using 



A nontrivial check of the calculation is then provided by the expected equality, for 
each diagram involving massless particles only, of the coefficient of the 2/e ultraviolet 
divergence and of the coefficient of the — log q 2 contribution, which would appear 
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as one single logarithmic term log(A 2 /<^) in the presence of an explicit ultraviolet 
cutoff A. 




2a) 



2b) 




2g) 

Fig 2. Additional one loop graviton exchange diagrams. 
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3 Results and Discussion 



By converting the expressions for the individual diagrams to coordinate space, one 
obtains the following results. One has from diagram la 



3 2 m 1 m 2 (m 1 + m 2 ) 2 m x m 2 
+ 7 u i h z u 3~' 



from diagram lb 
from diagram lc 

from diagram Id 
from diagram 2b 
and from diagram 2d 



3 2 m l m 2 {m l + m 2 ) 2 mim 2 
+ — Lt h Z It — . 



2 m 1 m 2 (m 1 + m 2 ) 2 ™i™2 
(j h o (j 



10 G 



2 mim 2 



16 2 ^1^2 
t — (j ; — • 



23 2 ™i™2 

+ "T U 3~ 

3 wr* 



(3.1) 

(3.2) 

(3.3) 
(3.4) 
(3.5) 
(3.6) 



From diagrams 2e and 2g one obtains the graviton and ghost vacuum polarization 
contribution 

(3.7) 



206 2 m i m 2 
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This last contribution was also computed in Ref. [12]. We have verified that the 
Slavnov- Taylor identity for the vacuum polarization H a f) 1 s{q), 



q^q v D l _ t \ a fj(q) Hafj-zsiq) D 1 s ua (q) — 



(3.8) 



is indeed satisfied to this order. In Ref. [4] the vacuum polarization was computed 
using a somewhat different expansion for the metric field, and a coordinate invariant 
expression for the one-loop counterterms was given in terms of operators quadratic 
in the curvature. 

Finally, diagram 2h represents the contribution to the vacuum polarization due 
to one massless scalar particle, 



1 2 m i m 2 



(3.9) 



20 irr- 

Its contribution to the vacuum polarization satisfies separately the Slavnov- Taylor 
identity, as one would expect from the covariant conservation law for the energy- 
momentum tensor associated with matter. Diagrams 2a and 2c do not give rise to 
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any classical relativistic or quantum correction, while diagram 2f vanishes identically 
in dimensional regularization. Diagrams 2b, 2d, 2e, 2g and 2h give only quantum 
mechanical corrections, involving closed graviton loops in all cases, except 2b. 
The sum of all contributions from diagrams la to 2g is therefore 

1 miraofrai + m 2 ) 122 mim 2 

+ - G 2{ y V + — — G 3.10 

2 r l 15 7rr J 

The contribution of n species of massless scalar particles to the vacuum polarization 
(arising from diagram 2h) changes the quantum correction to the potential to 



1 , , (i mi mo 

+ - 488 - 3n G 2 ^^, 3.11 
oO 7rr J 



which represents a relatively small modification to the result for pure gravity if n 
is small. Massless particles of higher spin will contribute additional terms to the 
vacuum polarization. 

When the appropriate powers of c and % are put back in, one obtains the following 
final answer for the corrected potential in pure gravity, valid to order G 2 

y (r) m _ G !2mU_ «('"■+ - ,3.12) 

r 1 Zc z r Lbirc^r^ 1 

As we alluded to previously, two very different length scales enter in the correction 
to the static Newtonian potential, namely the Schwarzschild radii of the heavy 
sources, 2Gm 8 /c 2 , and the Planck length (Gh/c 3 ) 1 / 2 . As a consequence there are two 
independent dimensionless parameters that appear in the correction term, involving 
the ratio of these two scales with respect to the distance r. Presumably the above 
calculation is meaningful only if these two length scales are much smaller than the 
distance r. 

Our calculations are similar in spirit to the work of Ref. [6]. There the starting 
point is also a calculation of the scattering amplitude in the limit of small momentum 
transfer. The potential is defined there as the non-relativistic limit of the one particle 
reducible graphs in the crossed channel, which represents therefore a subset of the 
graphs considered here. We should point out that the results we obtain here are 
in complete agreement with the expected classical relativistic correction, as derived 
for example from the expansion of the Schwarzschild metric [20]. The sign of the 
quantum correction is found to be the same as in Ref. [6], and the magnitude of the 
correction is comparable. The sign of the quantum correction we obtain indicate that 
gravitational interactions increase (slowly) with distance, which shows similarities 
with the evolution of the coupling constant in pure Yang-Mills theories, but differs in 
sign from the QED radiative corrections to the static Coulomb potential. This result 
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is also in agreement with the intuitive expectation that gravity couples universally 
to all forms of energy, and cannot be easily screened by quantum fluctuations. 

Recently the authors of Ref. [10] have computed the corrections to the static 
Newtonian potential following the method of Ref. [21], thus extending to the next 
order in G the calculation of Ref. [22]. In their work the radiative corrections to the 
potential are obtained by considering correlations between the action contributions 
from two heavy particle world lines, separated by a fixed geodesic distance. The 
results they obtain appear to correctly reproduce the classical relativistic correction, 
but arise from only a subset of two diagrams among the four which lead to the 
classical correction in Ref. [18]. In this last reference the ladder and crossed ladder 
diagrams give, using the same metric expansion, additional contributions which 
appear to be necessary in order to obtain the correct classical relativistic correction. 
These diagrams involve recoil of the massive particles, and have been neglected in 
the calculation of Ref. [10]. In our calculation we find that ladder and crossed ladder 
diagrams (la and lb), when carefully treated, contribute to the quantum correction. 
This probably explains why our results and the results of Ref. [10] differ in both 
sign and magnitude for the quantum correction. 

Let us conclude by mentioning that we have little to say about what might hap- 
pen to higher order in the perturbative expansion. In particular it is unclear if higher 
order corrections in G can still lead to finite corrections in the long distance limit, 
as was found above to lowest non-trivial order. Whether higher derivative terms 
or string theory is needed to control the ultraviolet divergences appearing at higher 
loops remains an open question [23]. As we pointed out before, another important 
omission in the present calculation is represented by the absence of a cosmological 
constant term. This term substantially modifies the propagation properties of gravi- 
tons already at tree level, and leads to new, momentum independent, vertices and 
Feynman rules for gravitons which were not considered here. 

Finally there is the issue of the non-perturbative definition of the Euclidean 
path integral for quantum gravity, which suffers from the problem of the unbounded 
fluctuations in the conformal mode, and for which an integration over complex 
conformal factors has been suggested, followed by an integration over conformal 
equivalence classes of metrics. In the framework of perturbation theory we did not 
have to deal with these difficult problems. 
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